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Abstract 

Gravitational lensing in metric theories of gravity is discussed. I introduce a generalized approx- 
imate metric element, inclusive of both post-post-Newtonian (ppN) contributions and gravito- 
magnetic field. Following Fermat's principle and standard hyphoteses, I derive the time delay 
function and deflection angle caused by an isolated mass distribution. Several astrophysical systems 
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for an observational detection. Actual measurements distinguish only marginally different metric 
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I. INTRODUCTION 



The principle of equivalence provides a firm foundation to any conceivable theory of 
gravity. On the other hand, the derivation of Einstein's field equations contains a strong 



element of guesswork. It is, there 
defined as theories such that ^, 



'ore, very interesting to test metric theories of gravity 
27| : i) space-time is a Lorentzian manifold; ii) the world 
lines of test bodies are geodesies; in) the equivalence principle in the medium strong form 
is satisfied. General relativity, Brans-Dicke theory and the Rosen bimetric theory satisfy 
these postulates. In these theories, the usual rules for the motion of particles and photons 
in a given metric still apply, but the metric may be different from that derived from the 
Einstein's field equations. The basic assumption of the existence of a dynamical space-time 
curvature, as opposed to a flat space-time of special relativity, still holds. 

Different metric theories can be compared with suitable tests. Bending and time delay 
of electromagnetic waves are two important phenomena predicted by theories of gravity. 
A comparison among general relativity and other viable theories of gravity can be led on 
the basis of higher-order effects. Intrinsic gravito-magnetism is such an effect. Mass-energy 
currents relative to other masses generate space-time curvature. In particular, the effect 
of the angular momentum of the deflector has been studied by several authors with very 



different approaches, see 



10 



11 



23, 



24, 



251 ] . The ppN corrections to the metric 



element have also to be considered. These second order contributions, in the cases of either 

Eravitational deflection of light or relativistic time delay, have been discussed, respectively, in 
IS l^Q and [3]. In particular, using a parametrized expansion of the metric element to 
second order in the Newtonian potential, Richter & Matzner Q, Q, Q, 20] have examined 
phenomena of light deflection by Sun and Jupiter. 

In this paper, I discuss deflection and time delay of light rays in the usual framework of 



gravitational lensing 



17 



221 ] . As shown in [4, 23], the standard assumptions of gravitational 



lensing, i.e. the weak field and slow motion approximation for the lens, allow us to consider 
higher-order approximation terms in the calculation of lensing quantities, so that a very 
general treatment of lensing phenomena can be performed. 

On the observational side, gravitational lensing is one of the more deeply investigated 
phenomena of gravitation and it is becoming a more and more important tool for exper- 
imental astrophysics and cosmology. The impressive development of technical capabilities 
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makes it possible to obtain observational evidences of peculiar metric theories of gravity in a 
next future and to test the degree of accuracy of the Einstein's field equations. Furthermore, 
observations of gravitational lensing phenomena could demonstrate the inertia-influencing 
effect of masses in motion. In fact, the gravito-magnetic field, predicted in 1896-1918, has 
not yet a firm experimental measurement. 

The paper is as follows. In Section m I introduce the generalized metric element in the 
weak field and slow motion approximation. Section HTT1 considers the Fermat's principle in 
stationary space-times. Sections IIVI and IS contain a derivation of, respectively, the time 
delay function and the deflection angle. In Section EH I consider light rays passing outside 
the lens; several astrophysical systems are discussed. Section IVm is devoted to a summary 
and to some final considerations. 



II. A GENERALIZED METRIC ELEMENT 



The comparison of metric theories of gravity with each other and with experiments can 
be performed in the slow- motion, weak field limit. Let us consider the approximate ppN 
metric element generated by an isolated mass distribution, with e nergy density p, negligible 
pressure, and energy current density f = pv\ d<c. As shown in |l5Ul8j]. the knowledge of 
light propagation to any given order requires knowledge of every component of the metric 
to the same order. We first consider the generalized central mass solution of general rela- 
tivity (Schwarzchild metric) in isotropic coordinates and expand it as power series in the 
small parameter to the ppN order. Then, we multiply the terms of this expansion by 
dimensionless parameters. This expression can be generalized to an arbitrary mass distribu- 
tion by replacing — with the standard Newtonian potential U, solution of AU = AnGp, 
U <C c 2 . It is U ~ 0(e 2 ), with e denoting the order of approximation. Finally, we intro- 
duce the non-diagonal components of the metric tensor generated by mass currents. We 
can write goi ~ — 4Vi, where Vi is the gravito-magnetic potential, solution of AVi = inGpVi, 
Vi ~ Uv ~ 0(e 3 ). The final expression for the approximate metric element is 



ds 2 



1 + 2^ + 2,^ 



c 2 dt 2 - 



U 3 (U 



rfx 2 



Sp — - — cat. 



Asymptotically, the metric reduces to the Minkowski's one. (3 and 7 are two standard 
coefficients of the post-Newtonian parametrized expansion of the metric tensor P, I^tJ. f3 



is related to nonlinearity of mass contribution to the metric; 7 measures space curvature 
produced by mass. In general relativity, it is (3 — 7 = 1; in the Brans-Dicke theory, (3 = 1 and 
7 = ii^. e and u are non standard parameters, e takes into account the ppN contribution to 
the metric [2J; /x quantifies the contribution to the space-time curvature of the mass-energy 
currents and measures the strength of the intrinsic gravito-magnetic field p. In general 
relativity, e = \i = 1. Additional terms in a parametrized expansion of the metric element 
can also be considered |l5, 

The approximated metric element just introduced cannot describe every conceivable met- 



rame effects, violations of 
. However, the metric in 



ric theory of gravity. In particular, it does not consider preferred 
conservation of four momentum and preferred location effects [2< 
Eq. (^) should be obeyed by most metric theories, with differences among them occurring 
only in the numerical coefficients. 

We will assume that during the time a light ray interacts with an isolated distribution 
of matter, the configuration of that matter does not change significantly. Then, the metric 
element can be considered as stationary with 

C/(x) ~ -G I -P^-d 3 x, (2) 



x — X 



and 



V(x) ~ —G [ ( PV)(X 'W (3) 
Js? I x - x I 

Our approximation is accurate enough to encompass all astrophysical tests that can be 
performed in the foreseeable future and makes it possible to describe gravitational lensing 
in quite general metric theories beyond the post-Newtonian order, without knowing any 
particular field equation. 

III. THE FERMAT'S PRINCIPLE IN STATIONARY SPACE-TIMES 



Within the assumptions of geometrical 



optics, Maxwell's equations have approximate 



solutions in generic curved space-times [22j. From the principle of equivalence, it turns 
out that light rays are null geodesies. They can be characterized by the Fermat's principle 
which states that a light ray from a source to an observer follows a trajectory, from among all 
kinematically possible paths, whose arrival time is stationary under first order variations of 
the path, 5r = 0. The Fermat's principle takes a version of particular interest in stationary 



space-times 



14 . l22j . We consider a metric whose components g a p are function of the 



spatial coordinates x % only (roman indeces label spatial coordinates). On a null curve, it is 

ds 2 = g a pdx a dx^ = 0; 



for the future-directed light ray, 



cdt = -^W + dh 



fi'oo 



(4) 



where dip = ( —g^ + 



90i90j \ J™i 



'/nil 



dx l dx J defines the spatial metric 



0. 



Let us consider an asymptotically flat space-time. The arrival time of a light ray, whose 
spatial projection is p, to an asymptotic observer is given by 



c J P \/9oo 9oo 



and the Fermat's principle states 



(5) 



(6) 



5 / ndlp = 0, 

where the spatial paths p are to be varied with fixed endpoints; n is an effective index of 



refraction defined as 



_ 9io i . 1 
n = e + — — . 

fi'oo v fi'oo 



(7) 



where e % = ^ is the unit tangent vector of a ray. This version of the Fermat's principle is 
formally identical with the classical one. 



IV. THE TIME DELAY FUNCTION 



Let us go, now, to apply the above results to our approximate metric element. The proper 
arc length is 



dh 



c 2 V4 2 \c 2 



dl 



eucb 



(8) 



where c?/ euc i = \/ §ijdx % dxi is the Euclidean arc length. Since Eqs. (00 IHJ), the total travel 
time reads 



n dh 



1 f I , M 

1- (1 + 7K + 



(9) 



/^ + 7(l-|) 



3 

+ 4« 



dl 



eucl • 



The time delay of the path p relative to the unlensed ray po is 



AT = i jjTn dip - J dZpj . (10) 
Equation (|TU|) can be expressed as a sum of geometrical and potential time delays 

AT = AT geom + AT pot . 

The geometrical time delay, 



cAT geom = / dip - / dip, (11) 

</p </p 

is due to the extra path length relative to the unperturbed ray p . 

The potential time delay AT pot is due to the retardation of the deflected ray caused by 
the gravitational field of the lens; it is defined as the difference between the total travel time 
and the integral of the line element along the deflected path, 

cAT pot = / ndlp- dip. (12) 
J p J p 

In usual lensing phenomena, the deflection angle of a light ray is very small and we can 



treat the lens as thin 
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Since the deflection occurs essentially in a small region near the 
deflector, the actual ray path can be approximated by combining its incoming and outgoing 
asymptotes. This trajectory is a piecewise-smooth null geodesies curve consisting of a null 
geodesic from the source to the deflector and one from the deflector to the observer. It is 
useful to employ the spatial orthogonal coordinates (£i,£2, 0> centred on the lens and such 
that the /-axis is along the incoming light ray direction e m . The three-dimensional position 
vector x can be expressed as x = £ + Ze in . The lens plane, corresponds to I = 0. 




With these assumptions, the geometrical time delay is 

(13) 



cAT D 



-geom - 2 



where D s is the distance from the observer to the source, D d is the distance from the 
observer to the deflector and is the distance from the deflector to the source; r\ is the 
bidimensional vector position of the source in the source plane. 

The potential time delay can be written as the sum of three terms, 

ATpot = ATg + AT™ + AT p p p t N (14) 



The first contribution contains the post-Newtonian correction to the time delay. To calculate 
the contribution of order G N to the lensing quantities, we need the path of the deflected 
light ray to the order G N_1 |8(. So, an integration along the line of sight, e in , is accurate 
enough to evaluate the pN contribution. This corresponds to the Born approximation. It is 

cans = -toi f udi 



pot — 2 1 cu eucl 

c Ao.s. 

~ -2(1 + 7)^/ ^'S(^')ln^^ + const (15) 
where £ is the projected surface mass density 

/ P(tl)dl 7 (16) 

J l.o.s. 

and is a scale-length in the lens plane. 

The second contribution to the time delay derives from the gravito-magnetic field, 

cAT GRM _ 4/i I Y e (17) 



C° 



Since the dragging of inertia 



the thin lens hypothesis 



frames induces a correction of order e , we can again assume 
12l l22| and integrate Eq. jm over the unperturbed ray ej n Is is 
easy to generalize the results in (2^; we get 

cAT GRM^ g/i G /• ^' S (£')( v .e ta )^(^)i n Jizii + const; (i 8 ) 

( v-e in)i.o.s. is the weighted average, along the line of sight, of the component of the velocity 
v orthogonal to the lens plane, 

Mai - Jwe.')y.o<g (19) 

The last term in Eq. is the ppN one. The first contribution to the ppN time delay 
derives from non-linear interaction of matter with space-time, represented in the approximate 
metric element, Eq. by terms which contain the square of the Newtonian potential. We 
have 



-.PPN(I) 



1 



cAT Pot = -4 



\-P + l(l- 1 -)+\e Ju 2 dl cucl (20) 

Since the ppN order is 0(G 2 ), the Born approximation breaks down, so that we have to 
consider the integration of the Newtonian potential on the deflected path, calculated at 
order 0(G), 

cAT ppN( 2 ) s _ (1+7) y udlmd _ ^ udl \ (21) 



Higher order corrections to the geometrical time delay also contribute to the ppN time 
delay. The difference in the distances of closest approach of the deflected and undeflected 
light rays represents a pN quantity: first-order corrections in the calculation of this difference 
are ~ 0(G 2 ) and induce a ppN contribution to the time delay. We will consider explicitly 
this third contribution in the Section IVI1 ATpQ t N cannot be expressed in terms of an integral 
of elementary functions and projected quantities, such as E. 

We remind that the time delay function is not an observable, but the time delay between 
two actual rays can be measured. 

V. THE DEFLECTION ANGLE 

To derive the deflection angle, we apply the Fermat's principle, Eq. (|5j). In the formalism 
of the previous sections, the Fermat's principle can be restated as: actual light rays, given 
the source position, are characterized by critical points of the total time delay, i.e. AT(£) 
is stationary with respect to variations of The lens equation is then obtained calculating 

V C AT(£) = 0; (22) 

we get 

T7 = ^- (23) 

ex. = — cV^ATp t is the deflection angle, i.e. the difference of the initial and final ray 
direction. Once again, the post-Newtonian and the gravito-magnetic contribution to the 
deflection angle have a simple expression. It is 

^ N te)-2(l + 7 )§ / d^Ete')-!^!-. (24) 

The equivalence principle, special relativity and Newtonian gravitational theory imply that 
a photon must feel the gravity field of massive body. They yield only the "1" part of the 
coefficient in front of Eq. ([24)1 . This accounts for the deflection of light relative to local 
straight lines. However, because of space curvature, local straight lines are bent relative to 
asymptotic straight lines. The contribution proportional to 7 in Eq. (J24|) is just the bending 
due to the space metric, described by the gu components of Eq. (JTJ); 7 measures, at the 
post-Newtonian order, the curvature generated by an isolated mass and varies from theory 
to theory. 
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The contribution of the gravito-magnetic field to the deflection angle is 



GRM 



(0 



G 

c 3 Jyt 2 



c^'£(£')(v-e in > Lo . s .(£' 



(25) 



The parameter \i tests intrinsic gravito-magnetism in conceivable metric theories of gravity 



23 



3- 



5]. For an analysis of the integral in Eq. (|25|). I refer to 

The ppN contribution to the deflection angle has not a simple form in terms of £(£). 



VI. THE NEARLY POINT-LIKE LENS 



Let us consider light rays passing outside a matter distribution of total mass M and radius 
R. We will refer to this situation as to the nearly point-like lens. If we place the origin of 
the coordinates at the centre of mass of the system, the dipole term vanishes identically j^] . 
The Newtonian potential, to order e 4 , reads 



where Uj 2 is the quadrupole term, 

Uj 2 



U=-G^r + Uj 2 , 
X 



jn M d2 3cos 2 ^-1 

<J2<~r-, — I7T-K ■ 



(26) 



(27) 



J 2 is the dimensionless coefficient of the second zonal harmonic and d is the angle between 
the symmetry axis and the field point x. Assuming the deflector to be symmetric about its 
angular momentum vector L, it is 



cos d 



Lx 

iLllxl 



The lensing quantities at the post-Newtonian order, Eqs. (fI3| 121 }) . reduce to 



and 



cAT pN 



a pN «) = 2(l+ 7 ) 



-2(l+,)^,n(i 



GM £ 



(28) 



(29) 



(30) 



c2 e - 

The first higher order correction derives from the quadrupole moment. Since Uj 2 is a 
higher-order term, we can integrate along the unperturbed path. It is 

3 (L lo J + L^i + L 2 i 2 f 1 



j U j2 dl ~ l -J 2 GMR 2 J 
= -J 2 GMR 2 



L 2 
2(i-L ± ) 2 



(p + e? /2 



di 



L 



l.o. s. 



1 



(31) 



L is the modulus of the angular momentum; Li. . s . is the projection along the line of sight; 
L\ and L 2 are the components on £1 and £2, respectively; Lj_ = (-^-, 0) is the projection 
of the versor of the angular momentum on the lens plane. The contribution to the time 
delay is, see also j^J, 

l-2(£-L ± ) 2 - (L 



cAT, 



■h 



(1+7) 



GM 



pot — V- ' U „2 



l.o. s. 



1 



(32) 



The predicted deflection due to the second zonal harmonic reads, see also 

2 



0, 



-2(1 + J2R 2 



1 - 4(£-L, 



L 



l.o. s. 



£ + 2(£-L ± )L ± }- 



1 



(33) 



Circular symmetry is broken. 

Let us now calculate the ppN correction. The first contribution comes from those terms 
in the approximate metric element which are quadratic in U. Since 



U 2 dl 



cucl ~ (GM) 



, , 00 1 „ (GM) 2 

' 1 dl = 7T- 



observer 



00 



it is 



e + 1 2 

3 



4 



(34) 



-/5 + 7(l-|) + 



/ GM\ 2 1 



(35) 



3 

2 '~ ' ' V" 2/ ' 4 
The second contribution to the ppN time delay derives from the integration of the Newtonian 
potential on the deflected path. Since the linearized metric to the pN order is spherically 
symmetric, the photon path lies in the plane through source, lens and observer. As well 
known, the deflected trajectory at the pN order reads j^, |3] 



A^(l + 7 )^i(Z + (Z 2 + a i/2 ) 



(36) 



where A£ is the distance of the photon from an axis, parallel to the optical axis, passing 
through On the deflected path, the modulus of the distance of the field point reads 

A£ 



We have 



Udl 



Udl 



l.o. s 



(Z 2 + £ 2 ) 1/2 



4 



(37) 



+ 



(Z 2 + £ 2 ) 3/2 Z 2 + £ 



dl 



-(1 + 7) 



(38) 
(39) 
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the second contribution to the ppN time delay reads 



cAT (2) 



(l+7) : 



( 



GMYir 



(40) 



V c 2 ; e' 

The last contribution to the ppN order comes from higher order corrections in the estimate 
of the length of the deflected path. It is 

1 fdAC 2 ~ 



dl 



dl 



1 + 2 



+ 



I 2 



(12 + £2)1/2 /2 + £2_ 



(41) 
(42) 



integrating from —D/2 to +D/2, D ^> £, we get 

2 



dl ~ D 



l + (l + 7)' 



/ GMV 1 



:(l + 7) : 



/GM\ 2 7T 



(43) 



The last term in the right hand side of Eq. (|4H|) does not depend on £) and represents a local 
effect; it contributes to the ppN time delay, 



cAT 



(3) 



ppN 



i(l+7)' 



fGM Y ;r 



(44) 



Then, adding the three contributes, we get 



cAT p p f = vr 



and 



a ppN (£) 



7T 



2(l+ 7 )-0 + -e 



2(l + 7 )-/?+-e 



V 



c 2 ; r 



GM 



e 2 ' 



(45) 



(46) 



For /3, 7 and e of order of the unity, nearly 80% of the bending come from the standard 
post-Newtonian parameters (3 and 7, the remaining 20% arises from the non-standard e 
coefficient. Our estimate agrees with the earlier results in [?], Isl, Q ■ 

Eqs. (1451 EH)) can be easily generalized. For a not specified mass distribution, we have 

2 



cAT r 



PPN 



pot 



and 



&pp n (£) 



2(l+ 7 )-/3 + -e 



2(l + 7 )-/? + -e 



l.o. s. 



l.o. s. 



2 I dl enc ] 



€1? 



eucl- 



(47) 



(48) 
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In order to compare the effect of dragging of inertial frames on the deflection angle with 
the ppN contribution, we use the results for a finite homogeneous sphere. For a deflector 
rotating about the £ 2 -axis with angular momentum L, it is, outside the lens radius 



24, 



25], 



AT GRM M) = _^G L COS^ 

c 6 f 



and 



AC J 

a™(e,e) = cos 20, (50) 

AC T 

a^{U) = sin 20, (51) 

where £ and # are the polar coordinates in the lens plane. The gravito-magnetic field 
breaks the circular symmetry. Both the ppN and the gravito-magnetic contributions to 
the deflection angle decrease as £ -2 . The sign of the gravito-magnetic correction varies on 
opposite sides of the lens, so that it can be separated experimentally from other terms. 

For an angular momentum directed along the ^-axis, the quadrupole contribution to the 
time delay, Eq. (|32|). reduces to 

a^/o / \ GM „ 2 cos20 
cAT p t ~ (l + 7 )_^J 2 i? 2 —^. (52) 

As can be seen comparing Eqs. (|49 [ 152 )1 . the angular dependence of the gravito-magnetic 
correction differs from that in Eq. ([52)1 . The sign of the quadrupole correction does not 
change in the equatorial plane. Furthermore, for a deflector with angular momentum with 
generic direction in th e sp ace, the gravito-magnetic effect depends only on the component 
of L in the lens plane |23[ , whereas L\ s enters the quadupole correction. 

The magnitudes of the different contributions to the deflection angle are considered by 
investigating real astrophysical systems acting as lenses. It is enough to use the values of 
the coefficient in general relativity, /3 = 7 = e — fi — 1. We will consider light rays in the 
equatorial plane (9 = 0). 

The post-Newtonian deflection angle for rays grazing the solar limb is 1.75 arcsec; a ppN 
is about 11 /iarcsec, where the contribution of the non-standard e coefficient is ~ 2 /iarcsec. 
Given the angular momentum of the Sun, L & ~ 1.6xl0 48 g cm 2 s _1 the gravito-magnetic 
correction is ~ 0.7 /iarcsec. Very Long Baseline Interferometry (VLBI) has improved the 
accuracy of the measurements of the deflection of radio waves by the Sun to the milli- 
arcsec level. This is not enough to measure the higher order ppN and gravito-magnetic 

12 



contributions, so that the parameters /3, e and // cannot be determined. However, strong 



constraints on 7 can be put. It is 7 = 1.000 ± 0.002 



2l| . an impressive confirmation of the 



prediction by general relativity. In Brans-Dicke theory, this measurement constrains the u 



parameter, u ~ 500. 

For an early type star, L = 10 2 L & {^f^ Q- For M = 1-4M , R = 1.1R Q and for a 
light ray grazing the star's limb, a pN ~ 2.23 arcsec, a ppN ~ 18 /xarcsec, a GRM ~ 0.10 milli- 
arcsec. The gravito-magnetic correction overwhelms the ppN one by an order of magnitude. 

The gravito-magnetic field becomes even more significant for a fast rotating white dwarf, 
where L ~ V0.2GM 3 R H- For M ~ M , R ~ lO" 2 ^, f ~ 6R, it is a pN ~ 29.2 
arcsec, a ppN ~ 0.1 milli-arcsec, and a GRM ~ 0.032 arcsec. In this case, the gravito-magnetic 
correction is quite important. 

Now, we want to apply our approximation to a galaxy acting as a lens. We take M = 
1O 12 M , R ~ 50 kpc and L ~ O.1M kpc 2 s _1 , as derived from numerical simulations 
It is a pN ~ 0.80 arcsec, a ppN ~ 2.2 //arcsec, a M ~ 0.16 milli-arcsec. The gravito-magnetic 
correction overwhelms the ppN one by two orders of magnitude. 



VII. SUMMARY 



The time delay function and the deflection angle for a single lens plane have been consid- 
ered in the framework of metric theories of gravity. I used an approximate metric element 
generated by an isolated mass distribution in the weak field regime and slow motion approx- 
imation, expanded to the ppN order, and with non-diagonal components which include the 
effects of gravity by currents of mass. 

Simple expressions for the ppN corrections to the lensing quantities have been derived. 
The nearly point-like lens has been used to consider several astrophysical systems. This 
approximation for the deflector is quite rough, but, some of the times, astrophysics can be 
tough. The gravito-magnetic correction and the ppN contribution to the deflection angle are 
of the same order for intermediate main sequence stars, like the Sun. For early type stars, 
white dwarfs and galaxies acting as lenses, the gravito-magnetic term overwhelms the ppN 
one. Differently from both pN and ppN contributions, which depend only on the impact 
parameter of the incident ray, the gravito-magnetic field induces a deflection which breaks the 
spherically symmetry, so that it offers the best perspectives for an observational detection. 
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However, a mass quadrupole moment can make the situation worse on the observational 
side. 

Ground based instrumentations, such as VLBI, or satellites, such as Hipparcos, can mea- 
sure deflection angles, respectively in the radio-wave regime and optical band, with accuracy 
of nearly milli-arcsec. Since the 7 parameter appears in the post-Newtonian expression of 
the lensing quantities, this accuracy put strong constraints on it. However, the other pa- 
rameters which enter the approximate metric element, that is the standard (3 term and the 
non-standard e and ji coefficients, need more accurate measurements. Lensing by fast ro- 
tating stars, such as white dwarfs, could give some hints on the dragging of inertial frames, 
whose strength is determined by the \x parameter. 

New generation space interferometric mission, such as SIM by NASA (scheduled for 
launch in 2009), should greatly improve the experimental accuracy, so that gravitational 
lensing could address, in the near future, two very interesting topics in gravitation: the 
detection of gravito- magnet ism and possible discrepancies of gravity from predictions of 
general relativity. 
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